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ABSTRACT 


A  statistical  model  for  obtaining  a  lower  100(1 -a)^ 
confidence  limit  on  system  reliability  was  developed  for  the 

J..,  ctment  o£   Navy  Special  projects  Office  in  Hay  19 o$.      This 

b  od  is  a  part  of  NAWEPS  OD  2930>.   Tie  accuracy  of  this 
lo  2l  is  examined  .„■  computer  simulation  for  various  failure 
time  distributions.   The  simulation  results  are  presented  and 
discussed.  **»>***>,.  «  - 
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CHAPTER  I 
INTRODUCTION 

In  the  last  few  years  reliability  has  become  of  primary 
concern  in  the  development  of  most  large  weapon  systems.   This 
growing  importance  of  reliability  has  produced  a  great  need 
for  methods  of  predicting  the  overall  reliability  of  these 
systems.   Because  of  the  size  and  complexity  of  present  day 
military  weapons  system,  it  is  necessary  for  major  contractors 
to  have  many  sub-contractors  located  throughout  the  country. 
This  fact  complicates  the  problem  of  formulating  a  method  for 
computing  lower  confidence  limits  on  system  reliability  where 
the  reliability  of  the  system  is  computed  as  the  product  of 
the  reliabilities  for  the  components. 

Recently  a  system  reliability  model  for  computing  a 
lower  confidence  limit  for  the  reliability  of  weapon  systems 
has  been  developed  for  the  Navy  and  is  contained  in  "Guide 
Manual  for  Reliability  Measurement  Program,'1  NAVJ3PS  OD  2930^, 
15  May  I965.   The  model  was  prepared  for  the  Department  of 
Navy  Special  projects  Office.   The  document  states  that  the 
model  can  be  utilized  by  all  contractors  for  subsystem  relia- 
bility measurement  and  by  the  Navy  for  weapon  system  relia- 
bility measurement.   The  model  permits  the  combination  of  test 
data  from  all  levels  (from  component  to  weapon  system)  and 
tyxoes  of  tests  into  meaningful  component,  equipment,  and 
subsystem  failure  rates  and  reliability  indices. 


An  underlying  assumption  to  this  statistical  model  is 
constant  failure  rate,  i.e.,  the  exponential  failure  law  is 
assumed  to  hold: 

a  =  e-^  (1) 


where 


H  =  reliability 

X   =  failure  rate  per  mission 

It  is  the  purpose  of  this  thesis  to  examine  the  accuracy 
of  the  statistical  model  contained  in  OD  29304.  The  investi- 
gation that  follows  will  examine  the  model  when 

(1)  the  failure  distribution  is  other  than  exponential; 
in  particular,  when  failure  times  have  a  Hayleigh 
distribution,  i.e.,  increasing  failure  rate. 
Appendix  1  contains  a  graph  of  the  Rayleigh  distri- 
bution and  its  failure  rate  function. 

(2)  the  failure  distribution  is  a  binomial  distribution. 

In  later  chapters,  it  will  be  explained  that  the 
lCO(l-a)  percentile  of  the  simulated  distribution  of  the 
lower  confidence  limit,  denoted  by  A,    t      and  the  true 
reliability  of  the  system,  R  ,  should  be  equal.   Therefore, 
the  difference  will  be  used  as  a  measure  of  the  accuracy  of 
the  model.   As  an  example  of  some  of  the  results,  consider 
the  system  with  k   components,  each  having  a  reliability  of 
.995.   Samples  of  sizes  of  50  and  500  components  were  tested 
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with  a  planned  test  time  of  2  mission  units.   The  true 
reliability  of  the  system,  R  ,  is  .9801.   The  80th  percentile 
points  of  the  distribution  of  the  lower  confidence  limit  for 
Raylei&h  distributed  failure  times  were  computed  to  be: 

a         nl       Al  _  a       |A1  -  a  "  Rs[ 

.20        50        .957  .023 

.20       500        .961  .019 

Jaen   the  failure  times  were  binomially  distributed  and  the 
test  time  was  one  mission  unit,  the  80th  percentile  points  of 
the  distribution  of  the  lower  confidence  limit  were  computed 
to  be : 


a 


ni       Al-a       Al-cc"R 


.20        50       .9683         .0120 
.20       500       .9800         .0001 

Chapter  II  contains  a  complete  explanation  of  the 
statistical  model  in  0D  2930^!-.   Chapter  III  discusses  the 
simulation  procedures  and  presents  the  simulation  results. 
It  also  discusses  the  two  methods  used  to  determine  trie 
accuracy  of  the  model.   Chapter  IV  discusses  and  summarizes 
the  results  and  states  the  conclusions. 


CHAPTER  II 

EXPLANATION  OF   STATISTICAL  MODEL  (1) 

If  a  system  consists  of  k  components  In  logical  series, 
the  true  system  reliability,  R  ,  may  be  expressed  as 

k 
R     =    TT  \  (2) 

s    1  =  1 


..here  R  is  the  true  reliability  of  the  lth  component, 
i 

It  is  desired  to  find  a  lower  100(1  -  a)%   confidence 

A 

limit  for  Rt.   That  is,  a  statistic  R^   ,  ,  such  that  what- 
ever the  actual  values  of  R  ,  i  =  1,2 , . . . ,k, 

pERs»i,L(a)]  =  1-a  (3) 

The  proposed  method  to  find  this  lower  100(1  -a)^ 

confidence  limit  for  R  as  outlined  in  OD  29304  will  be  ;re- 

s 

sented  below  to  include; 


(1)  Operating  assumptions. 

(2)  Estimation  of  failure  rates. 

(3)  Statistical  equations  for  confidence  limits  on 
failure  rates  and  reliability. 

Op  e  r  a  1 1  n^,  As  a  u  jip  t  i  o  n  s 

The  reliability  Measurement  system  described  In  OD  29304 
is  predicated  upon  the  general  reasonableness  of  the  follow- 
ing assumptions: 
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(1)  Constant  Failure  Rate  -  the  exponential  failure  law 
is  assumed  to  hold. 

(2)  Additivity  of  Stress  Effects  -  the  failure  rate  in- 
duced "by  two  simultaneously  acting  stresses  is 
equal  to  the  sum  of  the  failure  rates  due  to  the 
two  stresses  acting  sequentially.   This  assumption 
permits  adding  data  from  single  environment  tests 
together  to  simulate  mission  experience. 

(3)  Independence  of  Component  Failures  -  independence 
of  component  failures  is  assumed  because  components 
are  normally  tested  individually  by  type,  and  sub- 
system reliability  is  estimated  using  component  and 
other  applicable  test  results  (e.g. ,  equipment  and 
subsystem) . 

(^■)  Failure  Rate  Constancy  -  the  failure  rate  is  con- 
sidered a  function  of  only  the  stress  acting.  In 
other  words,  the  item  under  test  has  no  memory. 

Estimation  of  Failure  Rates 

Under  the  above  assumptions,  if  component  1  has  failure 
rate  A.-,,  and  component  2  has  failure  rate  x.? ,  then  the 
corresponding  failure  rate  for  the  simple  serial  subsystem 
is  : 

^■subsystem  =  xi    ''   kz  ^ 

Since  the  mode  of  testing  generally  employed  by  subsystem 

contractors  is  to  test  until  some  planned  test  time,  the 

A 

failure  rate  estimates,  \    ,  contain  bias.   Therefore,  an 

l 
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unbiased  failure  rate  estimate  has  been  developed  In  OD  2930^ 
and  is  given  as 


i   Sj. 


2ni 


2n  +  1 


(5) 


where 


iii  =  sample  size  for  the  ith  component;  i.e.,  the 
number  of  tests  conducted  on  component  i. 

S  =  the  sum  of  all  test  times  accumulated  on  the 
n  components  of  type  i. 

f  =  number  of  components  of  type  i  that  did  not 
complete  their  mission,  i.e.,  failure  time 
was  less  than  the  planned  test  time. 

Under  the  OD  2930^-  method,  failure  rate  estimates  for 
subsystem  reliability  can  be  calculated  from  component, 
component-environment,  or  component-environment-test  con- 
dition level  data.   That  is 

A  A 

^•component  ~  £  ^component-environment        (°) 
A  A 

^•equipment  ~  ^  ^component  *'  ' 

A  A  .q  . 

subsystem  '  L     equipment 

It  permits  the  combination  of  test  data  from  all  levels 
(from  component  to  weapon  system)  and  types  of  tests  into 
meaningful  component,  equipment,  and  subsystem  failure  rates. 
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Statistical  Equations 

All  statistical  estimates,  based  as  they  are  upon 
fragmentary  data,  are  subject  to  statistical  uncertainty. 
In  xeliability  work,  it  is  usually  of  interest  to  compute 
upper  limits  on  the  failure  rate  and  corresponding  lower 
limits  on  the  reliability.   The  following  statistical  equa- 
tions for  calculation  of  the  lower  confidence  limit  are 
based  on  normal  theory  and  have  been  corrected  to  compensate 
for  small  values  of  x.      Detailed  derivations  of  these 
formulas  are  contained  in  reference  1. 

The  method  proposed  is  to  put  n   items  of  component 
i  on  test  under  the  environmental  conditions  defined  in  the 
mission  and  let  each  operate  until  failure  or  the  planned 
test  time,  whichever  occurs  first.   All  time  is  measured  in 

mission  units.   Then  an  unbiased  estimate  of  failure  rate, 

A 

X.   ,  for  each  component  i  is  obtained  from  equation  (5). 

If  failures  have  occured  for  any  of  the  i  components ,  an 
upper  limit  on  failure  rate  can  be  calculated  as  follows: 


2* 


2A         ,       ,^A2 


\     r-.   2JL±  (3K)    C  +Nl4?i(|3K)    C  +    (gK)    G  ^ 


'u  2 


where 


A     k  A 

X   =   £  X,  (10) 

i=l 

A 


S-  StySj)  (U) 

A 

X 

K  =  (l-cc)th  percentile  of  the  normal  distribution 
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3  =  correction  factor.  Table  1  Shows  appropriate 
beta  values  for  the  80th  percent  confidence 
limit 

/hen  no  failures  have  occurred,  the  upper  limit  on  failure 
rate  is 


k 


Lh 


1  =  1 


w.iere 


k  =  number  of  components  in  the  simulations 

3  =  the  sum  of  all  test  times  accumulated  on  the 

i 

n   components  of  type1  i 
i 

Substitution  of  the  upper  limits  on  the  failure  rates 
previously  obtained  will  generate  corresponding  lower  con- 
fidence limits  on  reliability 5  thus: 

ft       =  e  Au  (13) 

s,L(a)  KXJ) 

To  determine  the  accuracy  of  the  above  model,  one  needs 

A 

to  look  at  the  characteristics  of  R      .   3y  def  iiiition, 
A  s,L(a) 

R    ,  ,  is  a  lower  100(1-  a )%   confidence  limit  for  a  ,  the 
s,L(cc)  s' 

system  reliability;  i.e.,  from  equation  (3) 

p[r  ,      <  a  ]  =  1  -  a 
L  s,L  1  ~  sJ 


This  says  that  R  is  always  the  (l-a)th  percentile  point 

S  A 

of  the  probability  distribution  of  R    .      .       Thus  if  we 

s  ,  L(oc ) 
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construct  the  distribution  of  R    ,  .  by  computer  simulation, 

s  , hoc  ; 

we  should  find  that  the  (l-a)th  percentile  point  of  our 
constructed  distribution  is  R  ,  if  equation  (3)  is  correct; 

A 

i. e.,  if  in  fact  R   ,    ,  is  a  true  100(1-  a )%   lower  confi- 

S  ,  L\CL  ) 

dence  limit  for  R  . 

s 

Thus  a  measure  for  the  accuracy  of  the  model  is 


|a_    -  r  I 

1  -  a    si 


where 


A     =  (l-a)th  percentile  Doint  of  the  distri^ 

bution  of  R 

s,L(a) 

If  equation  (3)  is  correct,  then  A,     should  equal  R  . 

1  -  a  s 


TABLE  1 
BETA  VALUES  FOR  80$  CONFIDENCE  LIMIT 


Number  of 

Beta 

Failures 

fi 

Values 

fi 

h 

f 
i 

*i 

0 

1.507 

10 

1.173 

20 

1.129 

1 

1.369 

11 

1.167 

21 

1.126 

2 

1.309 

12 

1.1.61 

22 

1.124 

3 

1.272 

13 

1.156 

23 

1.122 

J* 

1.246 

14 

1.153 

24 

1.119 

5 

1.227 

15 

1.146 

25 

1.117 

6 

1.212 

16 

1.144 

26 

1.115 

7 

1.200 

17 

1.139 

27 

1.113 

8 

1.190 

18 

1.135 

28 

1.111 

9 

1.180 

19 

1.132 

29 

1.110 

15 


CHAPTER  III 

SIMULATION  PROCEDURE  AND  RESULTS 

In  order  to  evaluate  the  accuracy  of  the  OD  29304  model, 

a  computer  program  was  written  and  used  to  simulate  the 

distribution  of  !EL   ,  *.   The  distribution  is  constructed  on 

s  ,n>i; 

a  computer  by  generating  500  values  of  R    .  .  for  a  Riven 

s ,L(a) 

set  of  parameter  values  k,  n  (i  =  1,2 , . . . ,k) ,  R  (i  =  1 ,2 , . . .  ,k) , 
Tq* (1  ■  1,2,. . • ,k),  and  a  given  failure  distribution;  i.e., 

A 

R    ,  ,  is  a  function  of 
s,L(a) 

k:   number  of  components  in  the  system 

n  :   number  of  items  of  component  i  tested 

R  :   reliability  of  component  i 
i 

T  . :   planned  test  time  for  component  i 

f,(t):   failure  distribution  for  component  i 
l 

.f\en   a  failure  distribution  for  component  i  is  specified, 

the  computer  program  generates  500  random  failure  time 

variates ,  T.  .  ( j  =  1,2  , . . .  ,500  ) ,  for  component  i.   If  the  gener- 

ated  time  variate  is  less  than  the  planned  test  time,  T   ,  a 

oi 

failure  is  "counted"  by  the  computer.   tfith  the  total  number 
of  failures,  f  ,  counted  for  a  set  of  500  generated  failure 

times,  the  computer  program  then  computes  trie  failure  rate 

A 

for  component  i,  \    ,  from  equation  (5).   The  above  procedure 

is  repeated  for  each  component  i,  and  the  corresponding 

A 

upper  limit  on  tne  system  failure  rate,  xu,  is  tnen  calculated 

from  equations  (9),  (10),  and  (11).   Vita  this  value  of  the 


lo 


system  failure  rate,  the  lower  confidence  limit  on  the 

A 

system  reliability,  R      ,  is  computed  from  equation  (13). 

s  »L(a) 

The  above  procedure  is  replicated,  giving  500  values  of 

A 

Ro  t t    \»   These  500  values  are  then  ordered  by  a  separate 
s  iHa ; 

subroutine,  and  the  computer  "picks"  the  (l-a)th  percentile 

of  these  500  values,  denoted  by  A,    .   This  value  is  then 

1  -  a 

compared  with  the  true  system  reliability  R  . 

s 

The  above  simulation  procedure  was  carried  out  for  tne 
following  two  methods: 

(1)  Rayleigh  distributed  failure  times 

(2)  Binomial  distributed  failure  times 

Rayleigh  Distributed  Failure  Times 

Since  the  Rayleigh  distribution  is  a  one  parameter 
distribution  (see  appendix  1),  it  was  decided  to  choose  the 
value  of  this  parameter  (a)  such  that  the  reliability  of  a 
component  for  one  mission  unit  is  .995;  i.e.» 

H1(l)  =  PC^j  >1]  =  .995  (1*0 

The  determination  of  the  value  of  this  parameter  Is  shown  in 
appendix  2.   The  resulting  failure  time  distribution  that 
will  be  assumed  for  each  component  i  in  the  simulation  is: 

f  (t,a)  =  ite   2a2     ,  a=10,  t>0       (15) 

1       or 

For  Rayleigh  distributed  failure  times,  the  accuracy  of  the 
model  was  examined  for  five  cases: 


17 


(1)  Case  1 

(2)  Case  2 

(3)  Case  3 


'oi 


ol 


'ol 


(k)      Case  4:   T 


oi 


(5)   Case  5:   TQl  = 


0.5  mission  units 

1.0  mission  units 

1.5  mission  units,  1  =  1,2 

1.35  mission  units,  1  ■  3,** 

2.0  mission  units 

5.0  mission  units 


Determining  the  Case  k   planned  test  time,  To1,  is  shown  in 
appendix  2.   For  each  case,  the  number  of  components  k  was 
equal  to  k,    the  reliability  HA  was  equal  to  .995»  and 
various  values  of  n^  were  used.   The  cases  are  numbered  and 

the  results  shown  in  Table  2,   The  mean  and  variance  of  the 

A 

distribution  of  R  r/  .  are  also  shown  in  Table  2.   The 

s ,L(a) 

quantity  labeled  TT  In  Table  2  is  expressed  by 


k 


TT=  £  n.  X 
1=1  *  1 


(16) 


where 


1  5  l  -  1  =  the  unreliability  of  component  i 

This  quantity,  TT,  is  a  function  of  the  amount  of  testing 
relative  to  the  unreliability  or  failure  rate.   The  results 
of  the  cases  in  Table  2  imply  that  the  model  is  more 
accurate  as  TT  increases. 

Binomial  Distributed  Failure  Times 

— ■— i  i   i— — — ——  i  ■   i— — "■  ■^—— — — —  — j— — 

A 

For  this  method,  the  failure  rate  for  component  i,  k*  » 
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is  given  as 


\  =*i  =  n;  <17> 


wnere 


f .  =  number  of  components  of  type  i  that  failed 
n.  =  number  of  items  of  component  i  tested 

Using  a  uniform  random  number  generator,  the  computer 
generates  a  uniform  random  value  betx^een  0  and  1  for  com- 
ponent i.   If  the  value  is  greater  than  R. ,  a  failure  is 

1 

"counted"  by  the  computer.   This  procedure  is  repeated  n 
times  for  each  R  and  for  each  component  i.   Thus  the  num- 
ber of  failures  counted  divided  by  n.  becomes  an  estimate 
of  the  failure  rate  for  component  i  as  given  by  equation 
(17)  above.   This  process  is  repeated  for  all  components, 
iith  these  estimates,  \j_,  the  computer  simulation  determines 

A.n     and  calculates  the  mean  and  variance  of  the  500 

1"a     A 

values  of  R    ,  .  as  previously  explained. 

S  ,  J_i  \  cc ) 

For  this  method,  the  accuracy  of  the  model  was  exam- 
ined for  three  cases: 


(1)  Case  1 

(2)  Case  2 

(3)  Case  3 

The  sum  of  the  failure  times,  S. ,  was  set  equal  to  to  n^;   i.e., 

ni 

S.  =  2  T,  .  =  n    ,  1  =  1,2 k  (18) 

1   i  =  1  1J     1 
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R 
i 

= 

.995 

R. 

l 

= 

.950 

% 

= 

.900 

the  number  of  components  k  was  equal  to  ^5  and  various 

values  of  n.  ,:ere  used.   The  cases  are  numbered  and  the 

1 

results  shown  in  Table  3.   Again,  the  quantity,  TT,  is 
also  listed  in  Table  3. 
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CHAPTER  IV 

SUMMARY  AND  CONCLUSIONS 

It  was  stated  previously  that  a  measure  of  the  accuracy 
of  the  model  was  given  by 


A.     -  R  I 
1  -  a    si 


where  A     is  the  (1-a)  percentile  point  of  the  distri- 

1  -  a 

bution  of  R_  r/~\,  and  R  is  the  true  reliability  of  the 

S  j  J_i  ^  CX.  /  s 

system.   For  the  Rayleigh  distributed  failure  times,  this 

measure  for  the  various  cases  is  shown  below  for  R  equal 

s 


o  .995. 

T  . 

01 

a 
.20 

50 

100 

500 

R 
s 

K        -   R  1 

1  1  -  a    s| 

0.5 

.995 
.995 
.995 

.0^3 
.012 
.008 

1.0 

.20 

100 

300 
500 

.995 
.995 
.995 

.010 

.003 
.003 

2.0 

.20 

50 

100 

500 

.995 
.995 
.995 

.023 
.019 
.019 

5.0 

.20 

50 

100 

50c 

.995 
.995 
.995 

.073 
.071 
.071 

The  above  results  bring  out  some  rather  interesting  points, 

If  the  planned  test  time,  T   ,  is  large,  the  model  gives 

oi 

pessimistic  results;  i.e.,  the  accuracy  is  not  very  good, 
regardless  of  the  sample  size  n. .   However,  if  the  planned 
test  time  is  small,  the  accuracy  of  the  model  is  very  good 
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if  a  large  sample  of  a  particular  component  is  tested.   In 
other  words,  for  small  planned  test  times,  the  accuracy 
of  the  model  is  very  good  when  the  quantity,  TT,  is  greater 
than  6 . 

For  the  Binomially  distributed  failure  times,  the 
accuracy  of  the  model  for  the  various  cases  is  given  below. 

a  Rq  n        IjL     -  R 

s  i        I  1  -  a    s 


.20        ,980  25  .043 

.980  50  .012 

.980  100  .010 

.980  500  .000 

.20        .950  25  .013 

.950  50  .008 

.950  100  .004 

.950  500  .001 

.20        .900  25  .013 

.900  50  ,004 

.900  100  .002 

.900  500  .001 


The  results  clearly  show  that  the  accuracy  of  the  model  is 
extremely  good,  especially  when  the  sample  size,  n. ,  is 
large. 

Obviously,  the  main  reason  for  the  results  In  the  case 
of  Rayleigh  distributed  failure  times  is  the  fact  that  the 
components  have  increasing  failure  rates.   The  underlying 
assumption  in  this  model  is  that  components  under  test  have 
a  constant  failure  rate.   In  practical  applications,  many 
components  follow  this  assumption  reasonably  well,  when 
infant  mortality  failures  are  deleted  and  wearout  failures 
are  not  anticipated  during  the  mission  duration.   However, 
if  a  new  component  is  being  tested,  a  constant  failure  rate 
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may  not  be  reasonable  at  all.   Therefore,  when  this  model 
is  being  used  to  obtain  a  lower  confidence  limit  on  the 
reliability  of  a  weapon  system,  and  it  is  known  that  a 
particular  component  has  an  increasing  failure  rate,  then 
the  points  brought  out  above  as  a  result  of  the  simulation 
should  be  taken  into  consideration. 
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APPENDIX  1 
HAYLEIGH  DISTRIBUTION 


PDF:   F(t)  ■  -^  t  e   ^a< 
a 


It2 
CDF:   F(t)  ■  1  -  e    * 


FH  (failure  rate)   z(t)  =  -—   (increasing  failure  rate) 

a2 


i 

) 

j06- 

•05 

J04- 

-03  J 

J02 

.01 

4   i 

a  =  10 


2        ^6         B      10      12    lb      16     IB      20     22     2U        t 
Probability  Density   Function 


E(tJ 


! 


.30 
25  - 

•20 
.15 
,10 

.05  1 


as  10 


%      $        6      8     10    "5     3f      2      18     20    22     2U~     t 
Failure   Rate   Function 
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APPENDIX  2 

Value  of  Parameter  for  Rayleigg  Distribution 

The  value  of  the  parameter,  a,  was  chosen  such  that 

R^l)  -  P[T>1]  -  .995  (19) 

where 

t2 


2  a2 
R.  (t)  -  1  -  F(t)  -  e  (20) 

Equating  equations  (19)  and  (20),  we  obtain 

1 
e  "a   =  .995   ,  t  =  1  (21) 


or 


a  -  "  2 (in  .995)  (22) 


hence 


a  =  10 


Planned  Test  Time  (T  . )  for  Case  k 
oi 

The  Case  k   planned  test  time  was  chosen  such  that  T 

oi 

satisfied   the    f ollovjing    equation: 


(23) 
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'daere 


2(t)  =  JL     i  (24) 

a 

T'lis  value  of  T  .  makes  t'ae  "average"  value  of  Z(t)  between 
01 

o  and  T  .  equal  2(1).   Setting  t  =  1  in  equation  (24)  and 


equating  (23)  ?.nd  (24)  gives 


or 


.ience 


oi 
\  t     It  s  -i  (25) 


T     2        2 

'o 


T   =20 

xoi     *L 


_1_ 

2 
a 
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APPENDIX  3 
COMPUTER  PROGRAM 


DIMENSION  PT(4,3).RSL(500),A(3).RSBAR(3),3RS(3)» 
3ETA(30),RS1(3),RSLBAR(3),RS(5C0).0(4),NM(4) 

READ(5,10)(NM(J),J=  1,4),  (BETA(h),li  =  1 ,30  ) ,  (  U  I ) ,  I  =  1  ^  ) 
10  F0RMAT(4i4/i6f5.3/14F5.3/4f6.3) 
TPR  =  URN(O) 
DO  999  K  =  1,4 
SS  a  NM(R) 
SRS(K)  =0,0 
A(K)  =  o5o 
RSLBAR(K)  =0.0 
SRSL(K)  =  0.0 
P2  =  NM(K) 
B  =  (2.0  *  F2)/(2;0  %  F2)  +  1.0) 

This  loop  generates  500  time  variates  for  each  component 

DO  900  M  =  1,500 

FS  =  0.0 

VU  =  0.0 

C  =  0.0 

VI  =  0,0 

RSL(M)  =0.0 

RS(M)  =  CO 

DO  800  I  =  1,4 

FF  =  0.0 

VL  =  0.0 

NL  =  NM(K) 

DO  700  J  =  1,NL 

TPR  =  URN(l) 

IF(TPR.GT.  2(I))GO  TO  605 

GO  TO  700 
605  FF  =  FF  +  1.0 
700  Continue 

VL  =  FF/F2 

VI  =  VI  +  VL 

FS  =  FS  4-  ^b" 

C  =  C  +  (VL/SS) 

VU  =  VU  +  (1.0/SS) 
800  CONTINUE 

This  statement  computes  a  value  of  R  =  e 

RS(M)  =  EXP(-Vl) 

RSBAR(K)  =  RSBAR(K)  +  RS(M) 

IF(FS  tGT.  0.0 )GO  TO  842 

BTA  =  BETA(L) 

VU  =  (VU»BTA)*(BTA»,1 772 41) 

Go  to  832 
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a 


8^2  C  =  C/Vl 

IP(PS  ,GE.  30.0)GO  TO  820 

LI  =  FS  +1.0 

BTA  =  BETA (LI) 

GO  TO  825 
820  BTA  =1.0 
825  B2  =  (BTA  *  .708964)  *  BTA 

D=  (('KO  *  V1)*(B2*C))  +  ((B2*B2)*(C*C)) 

D  =  S}RT(D) 

VU  -  (((2.0  «  VI)  +  (B2*C))  +  D)/(2.0)    A 

a  -X 

This  statement  computes  a  value  of  R  r  /„  \  =  e 

832  RSL(M)  =  EXP(-VU) 

RSLBAR(K)  =  RSLBAR(K)  +  RSL(M) 
900  CONTINUE 

The  following  statements  commute  the  mean  and  variance  of 

A        A 

R^  and  R 
s      s,L 

RSBAR(K)  =  R3BAR(K)/500.0 
RBAR  =  RSBAR(K) 
DO  9C5  M  =  1,500 

SRS(K)  =  SRS(K)  +  (RS(M)  -  RBAR)**2 
905  CONTINUE 

SRS(K)  =  S^RT(SRS(K)/500.0) 
RSLBAR(K)  =  RSLBAR(K)/5C0.0 
RLBAR  =  RSLBAR(K) 
DO   2^10  J  =  1,500 

SRSL(K)  =  SRSL(K)  +  (RSL(J)  -  RSBAR)**2 
2^10  CONTINUE 

SRSL(K)  =  S3RT(SRSL(K)/5C0.0) 

A 

The  following  statements  order  the  500  values  of  R  L  and 

picks  the  80th  percentile  of  the  500  values,  A-,   .' 

II  =  1 

DO  910  I  =  1,101 

TEMP  =  RSL(l) 

DO  915  J  =  1,500 

IF(TEMP  -  RSL(J))916,915,915 
915  TEMP  =  RSL(J) 

II  =  J 
915  CONTINUE 

RSL(II)  =  0.0 

A(K)  =  TEMP 
9.0  CONTINUE 

The  following  statements  print  out  the  results-  A,    ,  mean 

and  standard  deviation  of  Ro  ,.  ,  and  the  mean  and  sfaftdard 

deviation  of  R  .  ' 

s 
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JRITE(6,1010)  A(K),fiSLBAR(K),SfiSL(K)9RSBAH(K),SRS(K) 
1010  FORMAT (//6P13. 5) 

999  CONTINUE 
END 
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